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NASA Technical Note D-4727 

PRELIMINARY ANALYSIS OF A TITANIUM ALLOY HONEYCOMB 
SOLAR ABSORBER HAVING BLACKENED WALLS 

By William J. Bifano 
August 1968 

In calculating the apparent hemispherical emittance of a honeycomb -solar absorber, 
it was necessary to solve nonlinear integro-differential equations (eq. (9)) which describe 
the energy transport within the honeycomb. For TN D-4727, numerical solutions were 
obtained using a digital computer and a specially devised iterative technique which was 
formulated in reference 5. Subsequent investigation indicated that this technique was in- 
adequate for length -to -diameter ratios (L/d) of 3 or more. Hence, the results presented 
in TN D-4727 are incorrect, the error being progressively larger with increasing L/d. 
Consequently, a new computer program was devised which yielded correct solutions in- 
dependent of L/d. While the basic conclusions of the report remain unchanged, the 
corrected values of apparent hemispherical emittance, the critical performance param- 
eter for honeycomb solar absorbers, are somewhat higher. 

The attached corrected copies of figures 2 to 6, 8, and 9 and table HI should be in- 
serted in place of those presented in the report. 

In addition, equation (11) on page 6 should read 
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TABLE m. - CALCULATED VALUES OF APPARENT HEMISPHERICAL EMITTANCE 
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[Wall thickness, 0.003 in. (0.00761 cm).] 
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Figure 4. - Effect of thermal conductivity on apparent hemispherical emittance. 
Base temperature, 2060° R (1144 K); diameter, 0. 5 inch (0. 127 cm); cell thick- 
ness, 0.003 inch (0.00761 cm). 
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Figure 6. - Apparent hemispherical emittance as function of heat -conduction parameter. 



Figure 8. - Comparison of cell all temperature ratio for collimated and divergent 
solar input conditions for 0. . astronomical-unit distance from Sun. Refer- 
ence conditions: base temper ire, 2060 u R (1144 K); diameter, 0.5 inch 
(1. 27 cm); thermal conductivity, 0. 667 Btu per hour-inch-°R (0. 138 W/(cm)(K)) ; 
cell wall thickness, 0.003 inch (0.00761 cm). 
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ABSTRACT 

A titanium alloy hexagonal honeycomb with blackened walls is considered as an 
absorber of collimated solar energy. Circular cylindrical geometry is assumed for the 
analysis as an approximation of the hexagonal cell structure. The apparent hemi- 
spherical emittance of such an absorber positioned over a black surface is calculated. 
Results are presented for cell length -to -diameter ratios of from 1 to 7 and cell diam- 
eters of 0.25, 0.5, and 1.0 in. (0.635, 1.27, and 2. 54 cm, respectively). The col- 
limated incident solar flux is assumed sufficient to attain base surface temperatures of 
1860° and 2060° R (1033 and 1144 K). The corresponding weight per unit area of such 
an absorber is also calculated. 
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PRELIMINARY ANALYSIS OF A TITANIUM ALLOY HONEYCOMB 
SOLAR ABSORBER HAVING BLACKENED WALLS 
by William J. Bifano 
Lewis Research Center 

SUMMARY 

A titanium alloy hexagonal honeycomb with blackened walls was considered as an 
absorber of collimated solar energy. Circular cylindrical geometry was assumed as an 
approximation of the hexagonal cell structure for the analysis. Calculated values of the 
apparent hemispherical emittance were presented for base surface temperatures of 
1860° and 2060° R (1033 and 1144 K), cell length^to -diameter ratios from 1 to 7, and 
cell diameters of 0.25, 0.5, and 1.0 inch (0.635, 1.27, and 2.54 cm, respectively). 

The results indicated that for a titanium alloy honeycomb an apparent hemispherical 
emittance of the order of 0. 1 can be achieved with a cell diameter of 0. 5 inch (1. 27 cm) 
and a length-to -diameter ratio of about 6. The calculated weight per unit area of such 
an absorber fabricated from 0.006 inch (0.01523 cm) titanium foil is about 1.65 pounds 

9 

per square foot (0. 81 g/cm ), excluding the absorber base. 


INTRODUCTION 

Thermoelectric, thermionic, and Brayton- and Rankine -cycle turbine generators 
are examples of heat engines that could be used to provide auxiliary electrical power 
aboard future spacecraft. Such heat engines, using solar radiation as a source of 
energy, would require receivers capable of absorbing and retaining a large fraction of 
the incident solar flux for conversion to electrical energy. Because the efficiency of 
a heat engine is dependent on the temperature difference between the hot reservoir 
(receiver) and the cold reservoir (radiator), it is obviously desirable to operate at the 
highest possible receiver temperature consistent with thermal reradiation loss and 
material stability considerations. 

One method of utilizing solar energy is to collect the radiation with a large curved 
reflector and focus it into a cavity absorber. Another method of absorbing and retain- 


ing a large fraction of the incident solar radiation is to use a selective coating on the 
absorber surface. Such a coating would have a high absorptance for solar radiation 
(short wavelength region) and a low total emittance in the infrared region. A number of 
coatings having these properties have been tested (ref. 1); however, they have not yet 
been subjected to long-term testing at temperatures above about 1460° R (811 K). In 
addition, the performance of such absorber coatings is known to decrease with increasing 
temperature. 

A third technique for absorbing solar energy at high temperatures (above 1460° R 
or 811 K) is to place a honeycomb structure over the absorber surface, which would 
allow incident collimated solar energy to pass through to the absorber surface while 
limiting (attenuating) the reradiation of thermal energy from the absorber surface. 

Such a honeycomb structure would be fabricated from a material having a low thermal 
conductivity to minimize the conductive losses along the honeycomb wall (and subsequent 
radiation to space) and a high strength-to -weight ratio at elevated temperatures. 

Although the performance of honeycombs fabricated from fragile materials such as 
glass has been investigated for terrestrial use at temperatures below 1460° R (811 K) 
(refs. 2 and 3), no information was found in the literature for high -temperature honey- 
combs suitable for space application. This study therefore considers the performance 
of a titanium -foil honeycomb structure consisting of close -packed cells that are hexag- 
onal in cross section. The cell walls as well as the base surface on which the honey- 
comb rests are assumed black, that is, all incident energy is absorbed. The honey- 
comb structure is assumed to be oriented with respect to the collimated incident solar 
rays such that the cell walls receive no solar radiation. (The effect of the divergence of 
the solar rays is considered in the appendixes for one specific case at 0. 1 AU from the 
Sun.) The use of these assumptions results in a solar absorptance of 1; hence, only the 
radiant energy emitted from the honeycomb must be determined to characterize its 
performance. The ratio of the radiant flux per unit area emitted from a honeycomb at 
base temperature T b to that of a blackbody radiator at the same temperature is 
defined as the apparent hemispherical emittance 

Calculated values of apparent hemispherical emittance are presented for base 
surface temperatures of 1860° and 2060° R (1033 and 1144 K), cell length-to-diameter 
ratios of 1, 3, 5, and 7, cell diameters of 0.25, 0.5, and 1.0 inch (0.635, 1.27, and 
2. 54 cm, respectively), and a cell wall thickness of 0. 003 inch or 0. 00761 centimeters 
(total wall thickness of 0. 006 in. or 0. 01523 cm). The wall material thermal con- 
ductivity was fixed at 0. 667 Btu per hour-inch-°R (0. 1384 W/(cm)(K)) for most calcula- 
tions, although the effect of increasing the value to 1. 0 and 1. 33 Btu per hour-inch-°R 
(0. 208 and 0. 277 W/(cm)(K)) was also included for a representative case. 
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The weight -to-area ratio of a titanium (Ti)-6 aluminum (Al)-4 vanadium (V) alloy 
honeycomb structure, excluding the absorber base, is presented as a function of the cell 
length -to -diameter ratio for cell wall thicknesses of 0. 003 and 0. 001 inch (0. 00761 and 
0. 00254 cm). Additional information is also presented for the calculation of the apparent 
emittance for other suitable materials and/or cell dimensions, that is, for other com- 
binations of thermal conductivity, cell diameter, cell wall thickness, and absorber 
surface temperature. 


METHOD OF CALCULATION 
Description of Model 

For the determination of the radiant energy emitted from a cell of a hexagonal 
honeycomb (fig. 1(a)) for a given base temperature T^, the radiation emitted from both 
the base and the walls of the cell must be considered. The analysis was simplified by 
using a circular cylindrical cell with a diameter equal to the hexagonal cell minimum 
diameter (as shown in fig. 1(b)) to approximate, the hexagonal cell. The cylindrical cell 
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has a length L, a diameter d, a radius Tq, and a wall thickness t, which represents 
one-half the wall thickness of the close-packed hexagonal cells. Both the wall and base 
surfaces are assumed to be black (i.e. , all incident energy is absorbed), while the 
exposed edge of the cell wall is assumed to be perfectly reflecting. 


Procedure 

The analysis of such a circular cylindrical cell having nonisothermal walls is 
presented in reference 4 for the case of an external diffuse -radiation source and fixed 
uniform base temperatures. Since the case considered herein assumed a cylinder 
oriented with respect to perfectly collimated radiation, the analysis of reference 4 had 
to be modified. The modification reduced the problem to one in which the cylinder was 
considered independently of the external radiation; that is, the cylinder wall was 
assumed to receive no incident solar radiation. Again, as in reference 4, base surface 
temperatures were arbitrarily chosen by assuming that a heat balance existed between 
the base and the heat engine. The modified form of the analysis of reference 4 is pre- 
sented in the remainder of this section. 

From figure 1(b), a typical cylinder wall element located at x = x Q has a surface 
area cLAXq = (^d)dx Q . (Symbols are defined in appendix A. ) The energy balance on such 
an element is as follows: 

K‘L + K‘)c„„ d -° « 

If the temperature difference across the wall thickness is assumed to be negligible, the 
conduction term is 

where the thermal conductivity of the wall is assumed independent of temperature, and 

(j 7 d)t approximates the wall cross-sectional area. 

In the determination of the net radiation (dQ . j , all incident energy is assumed 

\ net /rad 

to be absorbed (i. e. , surfaces are black). The net radiation is therefore equal to the 
difference between the energy emitted by the element and the energy incident upon it. 

The energy emitted from the element, given by the Stefan-Boltzmann law, is 

eoT 4 (x Q )dA x (3) 


4 


p 


where e = 1 for a black surface. The incident radiation consists of contributions from 
both the cylindrical and base surfaces of the cavity. The radiation emitted by an element 
dA^. on the cylinder wall and absorbed by an element dA x is 


K) dF x-x 0 “A* 


(4) 


where dF x _ x ^ ^ diffuse -angle factor between an element dA x and an element 

dA • The reciprocity theorem for diffuse -angle factors states that dF dA = 
x o x-x o x 

dF „dA„ . Therefore, expression (4) becomes 
x o~ x x o 

oT 4 (x)dF _ x dA 

o o 


(5) 


The total radiation from the cylinder wall to the element dA is found by integration: 


dA x ^ 

°-4=o 


aT*(x)dF. 


x o“ x 


( 6 ) 


The radiation from the base surface to the element dA , found in a similar manner, is 

x o 


< ' T b F x -b dA x„ 
o o 


(7) 


Terms (6) and (7) are then subtracted from (3) to give the net radiation 

rL 


dQ 


neti 


rad 


o T 


*-L 


crT’(x)dF. 


_ - ctT^F^ . 
x -x b x -b 
o o 


dA. 


( 8 ) 


The energy balance equation (1) can now be evaluated. Substituting equations (2) and (8) 
into equation (1) and introducing dimensionless variables yield 


d 2 6 

dX 2 


= N. 


r L/d 
■/ ' 


e (X)dF Xo _ x - F v „ 


(9) 


where 
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d 



N c = 


^b 112 

kt 


The parameter N c is herein defined as the heat conduction parameter. Note that, for 
given values of the base temperature and cylinder dimensions, N c is solely dependent 
on the thermal conductivity of the wall material. 

Equation (9) is a nonlinear integro -differential equation from which the axial temper- 
ature distribution 6 along the cylinder wall can be calculated. A method similar to 
that used in reference 5 was used to solve this equation. The boundary conditions and 
angle factors are now defined. 

If the exposed edge of the cylinder wall is assumed to be perfectly reflecting, 
then at x = 0, dT/dx = 0. At x = L, the temperature T is set equal to T^. The 
boundary conditions are thus given by equations (10) as 


M = 0 at X = 0 
dX 

(10a) 

6 = 1 at X = — 

(10b) 


d 

The expressions for the angle factors, derived from angle -factor algebra, are 


r 


dF 


X o- X 


J, l x - X ol[ 2 ( X - X o) S+3 ] 

-x~ y rj? r 

€ x - X o ) 2 + *] 7 


dX 


(ID 


6 




( 12 ) 


Both equations (11) and (12) were obtained from reference 4. 

With a solution determined for the axial temperature distribution (eq. (9)), the 
amount of energy radiated from the cylindrical cavity can be calculated. From a given 
element dA x , the energy radiated out through the cavity opening is 

oT 4 (x)dA x F x _ 0 (13) 

The total amount of radiated flux emitted through the cavity opening from the cylinder 
wall is found by integration 

r L/d 

q w = 4(7T b / 0 4 (X)F x _ o dX (14) 

J 0 


where q,„ is the rate at which energy is radiated out of the cavity from the cylinder 
wall per unit area of the opening, and F Xq is the diffuse -angle factor between an 
element at X = X Q and a disk at X = 0. Adding the radiated flux emitted from the base 
surface to the flux from the wall gives the total rate at which energy is radiated out of 
the cavity per unit area of opening q^ 



' f L/i 4 

4 / Ax)F x . 0 dX + F b _ 0 

J o 


(15) 


where F^_ Q is the diffuse-angle factor between a disk at X = L/d (i. e., T = T^) and 
X = 0. The expressions for the angle factors are given as 


1 + 2X C 


X-o 


(x 2 + l) 


\ 1/2 


X (ref. 4) 


r b -o= 1 + 2 ''- 



+ l 


(ref. 6) 


(16) 


(17) 


7 



TABLE I. - VALUES OF HEAT -CONDUCTION PARAMETER FOR WHICH 
APPARENT HEMISPHERICAL EMITTANCE IS CALCULATED AT 
LENGTH-TO-DIAMETER RATIOS OF 1, 3, 5 AND 7. 


[Thickness, t, 0.003 in. (0.00761 cm).] 


Base temper- 
ature, T b 

Thermal conductivity, k 

Diameter of cylindrical cell, d, in. (cm) 

Btu/(hr)(in.)(°H) 

W/(cm)(K) 

0.25 (0.635) 

0.5 (1.27) 

1.0 (2.54) 

°R 

K 

Heat -conduction parameter, N c 

1860 

1033 

0.667 

0.1384 

2.413 

9.652 

38. 609 

2060 

1144 

.667 

.1384 

3.278 

13.112 

52. 45 

2060 

1144 

1.000 

.208 


a 8. 742 


2060 

1144 

1.333 

.277 


a 6. 556 



Evaluated for L/d values of 3, 5, and 7 only. 


Thus, , for given values of N c , L/d, and T b , it is possible to calculate the apparent 
hemispherical emittance e of the cavity, where 

a 

(18) 

With the use of a digital computer, numerical solutions of equations (9) and (18) were 
obtained for the values of N c and L/d shown in table I. Although the indicated values 
of the conduction parameter N c were determined for preselected values of T^, d, k, 
and t given in table I, the results are equally valid for other combinations of T b , L, 
d, k, and t yielding the given values of N c and L/d. 

RESULTS AND DISCUSSION 
Apparent Hemispherical Emittance 

The effect of varying the cell length-to -diameter ratio L/d on the axial temper- 
ature distribution is illustrated in figure 2, where the temperature ratio T/T^ is pre- 
sented as a function of the distance from the opening x/d for the following reference 
case: N c = 13. 112 (T b = 2060° R or 1144 K); d = 0. 5 inch (1. 27 cm); k = 0. 667 Btu 
per hour-inch-°R (0. 1384 W/(cm)(K)); and t = 0.003 inch (0.0076 cm)). Note that for 
an L/d of 5, for example, the wall temperature at the cell opening is only about half 


a 


q. r A , 

= _L = 4 / e 4 (x)F x _ 0 dx + F b . 0 

oT? •'O 
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Figure 2. - Temperature ratio as function of distance from cell 
opening. Base temperature, 2060° R (1144 K); diameter, 

0.5 inch (1.27 cm); thermal conductivity, 0.667 Btu per 
hour-inch-°R (0. 1384 W/(cm)(K)) ; cell wall thickness, 

0.003 inch (0.0076 'em). 



Figure 3. - Apparent hemispherical emittance as function of 
cell length -to-diameter ratio. Base temperature, 2060° R 
(1144 K) ; diameter, 0.5 inch (1.27 cm); thermal conduc- 
tivity, 0.667 Btu per hour-inch-°R (0. 1384 W/(cm)(K)) ; 
cell wall thickness, 0.003 inch (0.00761 cm). 


the base temperature. The effect of this temperature decrease in the region near the 

cell opening is to lower the apparent hemispherical emittance ~e , as shown in figure 3. 

At an L/d of 5, e a is 0. 13, and increasing L/d from 5 to 7 results in a decrease in 

e„ to 0. 065. 
a 

The thermal conductivity k of the Ti-6A1-4V alloy, obtained from reference 7, is 
presented in table n for various temperatures. Since the analysis assumed k to be 
independent of temperature, a reference value of 0. 667 Btu per hour-inch-°R 
(0. 1384 W/(cm)(K)) was used initially. Also, was evaluated for values of k of 
1. 00 and 1. 33 Btu per hour-inch-°R (0. 208 and 0. 277 W/(cm)(K)) and compared with 
the reference condition in figure 4 to assess the effect of variations in the wall material 
thermal conductivity. For an L/d of 5, for example, "c increases from 0. 130 to 
0. 155 when k is doubled from 0. 667 to 1. 333 Btu per hour-inch-°R (0. 1384 to 
0. 277 W/(cm)(K)). This effect of increasing conductivity might be offset by halving the 
wall thickness. 
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Cell length-to'diameter ratio, L/d 


Figured - Effect of thermal conductivity on apparent 
hemispherical emittance. Base temperature, 

2060° R (1144 Kh diameter, 0.5 inch (1.27 cm); 
cell thickness, 0.003 inch (0.00761 cm). 


TABLE H. - THERMAL CONDUCTIVITY OF SHEET 
TITANIUM - 6 -ALUMINUM -4 VANADIUM ALLOY 
[Data from ref. 7. ] 


Temperature 

Thermal conductivity, k 

°R 

K 

Btu/(hr) (in. )(°R) 

W/(cm)(K) 

672 

373 

0.417 

0.087 

860 

477 

.516 

.107 

1260 

700 

.716 

.149 

1660 

922 

.916 

.190 





p- 



0 1234567 0123 -4 567 


Cell length -to-dia meter ratio, L/d 

(a) Base temperature, 1860° R (1033 K). (b) Base temperature, 2060° R (1144 K). 

Figure 5. - Effect of cell diameter on apparent hemispherical emittance. Thermal conductivity, 0.667 Btu per hour-inch-°R (0. 1384 
W/(cm)(K)) ; cell wall thickness, 0. 003 inch (0. 00761 cm). 

The apparent hemispherical emittance ^ is presented in figures 5(a) and (b) for 
base temperatures of 1860° and 2060° R (1033 and 1144 K). The decrease in e a 
with increasing cell diameter d for a given L/d results from the attendant increase in 
cell length and thus from the reduced heat conduction. For cell diameters in the range 
of 0. 5 to 1. 0 inch (1. 27 to 2. 54 cm), ~e 's of the order of 0 . 1 result for L/d's of 6 or 

a 

more. In comparison, for an ideal selective solar absorber coating, as defined in 
reference 8, the theoretical hemispherical emittance is 0.03, and solar absorptance 
is 0. 87 for a surface temperature of 2060° R (1144 K) and a solar concentration of 10 
(ref. 8). However, such coatings for use at this temperature have not been developed. 

Values of the apparent hemispherical emittance of a cell, evaluated for a range of 
N c of 0 to 60, are presented in figure 6. Any combination of parameters T b , d, k, 
and t resulting in N c ’s in this range can be used to determine e a for a given L/d, 
as will be illustrated in the section Weight -to -Area Ratio of Honeycomb. 
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TABLE m. - CALCULATED VALUES OF APPARENT HEMISPHERICAL EMITTANCE 


[Wall thickness, 0.003 in. (0.00761 cm).] 


Diameter 

Cell length- 






Base temperature, T 

b 



of 

to -diameter 



— 








cylindrical 

ratio, 

1860 u R (1033 K) 




2060° R (1144 K) 



cell, 

d 

L/d 

Wall material thermal conductivity, k 

in. 

cm 


0.667 Btu/ (hr) (in. )(°R) 

0.667 Btu/ (hr) (in. ) (°R) 

1.000 Btu/(hr)(in.)(°R) 

1.333 Btu/ (hr) (in. ) (°R) 




(0.1384 W/(cm)(K)) 

(0. 1384 W/(cm) (K)) 

(0.208 W/(cm) (K)) 

(0.277 W, 

/(cm) (K)) I 




Heat- 

Apparent 

Heat- 

Apparent 

Heat- 

Apparent 

Heat 

Apparent 




conduction 

hemi- 

conduction 

hemi- 

conduction 

hemi- 

conduction 

hemi- 




parameter, 

spherical 

parameter, 

spherical 

parameter, 

spherical 

parameter, 

spherical 




N , 


emittance, 

N c 

emittance, 

N c 

emittance, 

N C 

emittance, 






? a 



e a 


? a 



0.25 

0.635 

1 

2.413 

0.715 

3. 

278 

0.684 







3 



.372 



.344 







5 



.212 



.189 







7 



.125 

1 


.107 





0.50 

1.27 

1 

9.652 

0.599 

13.112 

0.583 







3 



.279 



.268 

8.742 

0.283 

6.556 

0.297 



5 



.139 



.131 

8.742 

.142 

6.556 

.153 



7 



.070 


f 

.064 

8.742 

.073 

6.556 

.080 

1.0 

2.54 

1 

38. 609 

0.545 

52. 

45 

0.539 







3 



.245 



.241 







5 



.113 



.110 







7 

1 


.052 

\ 


.050 









Figure 6. - Apparent hemispherical emittance as function of heat -conduction parameter. 


Calculated values of e & for the specified conditions of T^, d, k, and t are 
presented in table HI. 


Absorber Efficiency 

In this analysis, the apparent hemispherical emittance is used to characterize the 
performance of a honeycomb solar absorber. Generally, however, the determination 
of the amount of usable solar energy made available by such an absorber is of more 
interest. The efficiency of the absorber tj is taken as'the ratio of the usable solar 
heat per unit area to that incident so that the absorber performance can be determined 
as a function of the distance from the Sun as follows: 

443 s' 2 - € 0T|| 
r\ = * — — 

443 s-2 


where 


o=l. 73xlO~ 9 Btu per hour-square foot-°R 4 
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TABLE IV. - ABSORBER EFFICIENCY FOR APPARENT 


HEMISPHERICAL EMTTTANCE OF 0. 1 AND BASE 
TEMPERATURE OF 1860° R (1033 K) 


Distance from Sun, 
s, AU 

Absorber 

efficiency, 

V, 

percent 

Usable solar heat per unit area 

Btu/(hr)(ft 2 ) 

W/m 2 

1.000 

0 

0 

0 

.465 

0 

0 

0 

.400 

.26 

720 

2 280 

.300 

.583 

2 870 

9 070 

.200 

.816 

9 050 

28 600 

.100 

.955 

1 42 250 

133 800 


or 


l 1 = 


0. 14 s" 2 - €„ot 5 
a b 


0. 14 s 


-2 


where 


-12 4 

o = 5.67x10 watt per square centimeter - K . 

Note that s is the distance from the Sun in astronomical units for both expressions, 
and a solar absorptance of 1 has been assumed. For the case of "e a = 0. 1 and 
= 1860° R (1033 K), the absorber efficiency as a function of the distance from the 
Sun is given in table IV. Note that the absorber efficiency is 0 from 1.0 to 0.465 astro- 
nomical unit, whereas from 0.3 to 0. 1 astronomical unit, it increases sharply from 
0. 583 to 0. 955. Hence, such a honeycomb absorber operating under the given conditions 
might be useful in the vicinity of Mercury (0. 388 AU) or at distances closer to the Sun. 


Weight-to-Area Ratio of Honeycomb 

The weight per unit area of a hexagonal honeycomb structure, excluding the ab- 
sorber base, is shown in figure 7 for Ti -6A1-4V alloy wall thicknesses t of 0. 001 and 
0.003 inch (0. 00254 and 0. 00761 cm). The physical density of the Ti - 6A1-4V alloy is 
0. 16 pound per cubic inch (4. 47 g/cu cm) (ref. 7). For an L/d of 6, for example, 
the weight -to -area ratio is about 1. 65 pounds per square foot (0. 81 g/cm 2 ) for 
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2.2 


LOi — 2.0 



0 1 2 3 4 5 6 7 

Cell length-to-diameter ratio, L7d 


Figure 7. - Hexagonal honeycomb solar absorber specific weight 
as function of cell length-to-diameter ratio for 6AI-4V titanium 
alloy with density of 0. 16 pound per cubic inch (4. 47 g/cu cm). 


t = 0.003 inch (0.00761 cm). Reducing t to 0.001 inch (0.00254 cm) at the same L/d 

O 

results in a proportionate decrease in weight to 0. 55 pound per square foot (0. 27 g/cm ). 

Titanium (6A1-4V) alloy honeycombs with a total wall thickness 2t as small as 
0.002 inch (0.00508 cm) have actually been fabricated, as reported in reference 9. 

Since the ?_ values were calculated for N„ values corresponding to a wall thickness 
t of 0. 003 inch (0. 00761 cm), figure 6 must be used to determine the effect on e. of 
reducing the wall thickness (changing N c ). For example, at an L/d of 7 and an N c 
of 3. 278 (when = 2060° R or 1144 K; d = 0. 25 in. or 0. 635 cm; k = 0. 667 Btu/ 
(hr)(in. )(°R) or 0. 1324 W/(cm)(K); and t = 0.003 in. or 0.0761 cm), the apparent 
hemispherical emittance e„ is 0. 106 in figure 5(b). Reducing the wall thickness to 
0.001 inch (0.00254 cm), that is, increasing N c to 9. 834, results in a decrease in 
€„ to 0. 07 as shown in figure 6. 

a 
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CONCLUDING REMARKS 


The analysis of the absorber performance presented in this report assumed a honey- 
comb structure perfectly oriented with respect to collimated solar radiation. Hence 
the effects of misorientation and the divergence of solar radiation should be considered 
for practical applications. Furthermore, a zero temperature gradient at the cell open- 
ing was used as a boundary condition by assuming that the exposed edge of the cylinder 
wall was perfectly reflecting. However, for certain cases, such as a distance of 

0. 2 astronomical unit or less from the Sun, the incident solar radiation absorbed at the 
cell edge would result in emittance values higher than those calculated. For example, if 
only the effect of the divergent radiation on the cell walls at 0. 1 astronomical unit from 
the Sun were considered for a length -to -diameter ratio of 5. 75 at reference conditions, 
the apparent hemispherical emittance e a would increase to 0. 19 as compared with an 
e_ of 0. 1 for collimated solar radiation (see calculations in appendix C). 

CL 

SUMMARY OF RESULTS 

A preliminary analysis of a titanium honeycomb solar absorber, with blackened 
walls positioned over a black surface, was performed for base temperatures of 1860° 
and 2060° R (1033 and 1144 K). Perfect orientation and collimated solar radiation were 
assumed. The study yielded the following results: 

1. With the incident solar flux assumed sufficient to obtain absorber surface tem- 
peratures of 1860° R (1033 K) or above, apparent hemispherical emittances of less 
than 0. 1 were indicated for values of length-to -diameter ratios greater than 6. 

2. A hexagonal close-packed cell structure was assumed for the honeycomb absorber 
fabricated from titanium - 6-aluminum -4-vanadium alloy foil 0. 006-inch (0. 0152 cm) 
thick. The weight-to-area ratio of such an absorber with a length-to-diameter ratio 

p 

of 6, excluding the absorber base, is about 1.65 pounds per square foot (0. 81 g/cm ). 
With the use of a foil thickness (2t) of 0. 002 inch (0. 00508 cm), from which honey- 
combs have been fabricated, the weight can be reduced to 0. 55 pound per square foot 
(0.27 g/cm 2 ). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 23, 1968, 

120-27-06-06-22. 


16 



APPENDIX A 


SYMBOLS 


A 

surface area 

0 

angle 

d 

diameter of cylindrical cell 

e 

hemispherical emittance 

F 

diffuse -angle factor 

t 1 

absorber efficiency 

h 

height 

e 

dimensionless temperature 

k 

thermal conductivity 


T/T b 

L 

depth of cell 

p 

variable radial coordinate 

N c 

heat-conduction parameter 

a 

Stefan-Boltzmann constant 

Q 

total heat-transfer rate 

01 

solid angle 

q 

heat-transfer rate per unit area 

Subscripts: 

r 

radius 

a 

apparent 

r O 

radius of cylindrical cell 

b 

base surface 

s 

solar radiation 

cond 

conduction 

s 

distance coordinate 

net 

net 

T 

absolute temperature 

0 

opening of cavity 

t 

cell wall thickness 

rad 

radiation 

X 

dimensionless axial coordinate 

s 

solar divergent rays 

X 

axial coordinate 

tot 

total 

a 

angle 

X 

axial coordinate 



w 

wall 
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APPENDIX B 


ESTIMATE OF SOLAR INPUT TO WALLS OF BLACK HONEYCOMB 
SOLAR ABSORBER AT 0. 1 ASTRONOMICAL UNIT FROM SUN 

Herein, the apparent hemispherical emittance e of a circular cylindrical cell with 
blackened walls was calculated with the assumption of perfectly collimated solar input 
to simplify the analysis. For the analysis, only radiant energy contributions from both 
the base and wall surfaces of the cell were needed to perform a heat balance on an element 
of the cylinder wall. Since the Sun is not a point source, however, the walls of a per- 
fectly oriented honeycomb absorber will receive a contribution from the divergent solar 
rays. This contribution will increase as the honeycomb distance from the Sun decreases; 
that is, as the solar rays become more divergent. The following approach is used to 
estimate the effect on ■?. of this additional heat input to the walls for a distance of 

a 

Oo 1 astronomical unit from the Sun: 



Assume the Sun to be a disk with radius r, area A 2 , and temperature T g of 
10 440° R (5800 K). The solar input to dA^ (see sketch) is determined by using the 
following approach: 


2,2 2 
s = h + p 

(Bl) 

dA 2 = p dp da 

(B2) 

cos p, -f> cos “ 
s 

(B3) 

1 _ 

cos /3 2 = — 
s 

(B4) 
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The heat a Q radiated per unit time from surface dA 2 within the solid angle under 


which dAj is seen from dA 2 is 


2 < 

d Q = — — cos /3g d^ 2 dA 2 
7 r 


(B5) 


where T g is the temperature of the Sun, 10 440° R (5800 K), and do> 2 is the solid angle 
under which dAj is seen from dAg. The expression 


dw 2 = 


dAj cos /3j 


(B6) 


gives the solid angle. Thus, 


<> crT cos /3, cos 
d^Q = — ? dAj dA 2 


(B7) 


(The heat radiated from dAj to dA 2 is negligible in this case.) Thus, 

/ 

'Ac 


dQ = 


ot; 


/ f cos 

\ 


cos /3r 


dA 2 dAi 


(B8) 


or 
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Evaluating equation (B9), where r is the radius of the Sun, 4. 33x10^ miles 
(6. 96xl0 10 cm), and h is the distance from the Sun, 0. 1 astronomical unit or 
9.3X10 6 miles (1.496X10 12 cm), yields 

-^§-= 1.99X10" 5 art 

/4A » 


= q o = 2. 8 Btu/ (hr) (in. 2 ) (0. 127 W/cm 2 ) 

A A S 
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APPENDIX C 

CALCULATED EMITTANCE OF BLACK HONEYCOMB SOLAR 
ABSORBER AT 0. 1 ASTRONOMICAL UNIT FROM SUN 

The results of appendix B are used to determine the effect on e 1 of an additional 
energy input to the walls of a honeycomb absorber at 0. 1 astronomical unit from the 

p 

Sun. The additional energy input is 2. 8 Btu per hour - square inch (0. 127 W/cm ). 


t — 

i— 


X T 

1 

dA* 

*0 

x o 

1 

L 

1 



1 


7 / 77777 / 


Consider a cavity with nonisothermal walls at 0. 1 astronomical unit from the Sun. At 

equilibrium, an energy balance on a typical element of the cavity wall located at x Q 

and having an area dA = (jrd)dx is 

x o ° 


K4o„d + K e ‘L =0 


(ci) 


The temperature variation across the wall thickness t is assumed negligible 


2 

(dQ A = -k( 7 rd)t dx where k / k(T) 
V ne Vcond ^2 0 


d^T 

= - kt — dA 


dx' 


The radiation emitted by the element is 


2 x o 


(C2) 


aT(Xo) dA x 


(C3) 
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The incident radiation comes from three sources: 

(1) Reradiation from the cylindrical surface of the body 

(2) Reradiation from the base surface of the cavity 

(3) Solar input due to divergent rays from the Sun 


Reradiation from Cylindrical Surface of Body 


The radiation leaving an element dA x on the cylinder wall and arriving at dA x is 


aT (x) dF x-x„ ^x 


(C4) 


where dF is the diffuse -angle factor between an element 
x-x o 

The reciprocity theorem for diffuse -angle factors states that 


dA„ and an element dA 
x x 0 


dF x-x. dA x = dF : 


x„-x dA x„ 
o o 


Therefore, egression (C4) becomes 




(C5) 


The total contribution from the entire cylindrical wall to the incident energy at x Q is 




-x 


(C6) 


Reradiation from Base Surface of Cavity 


The radiation from the base surface to dA is analogously 

x o 

<\-^ 0 


or 




Trdr 


(C 7) 


o 4 
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Solar Input Due to Divergent Rays From Sun 


The solar input to the wall from the divergent rays of the Sun at 0. 1 astronomical 

o 

unit is estimated to be 2. 8 Btu per hour - square inch (0. 127 W/cm ) = q . Thus, the 


total contribution due to the solar divergent radiation to dA is 

x o 

dA x„ 


(C 8 ) 


The net radiation from an element at x Q is then obtained by subtraction of expres- 
sions ( 6 ) to ( 8 ) from (3). Thus, 


dQ 


net 


rad 


<o> 


/ L 


"(x) dF x -x - orT b F x -b - «s 
x =0 0 0 


dA.. 


(C9) 


Substituting equations (2) and (9) into (1) yields 

°T 4 (Xo> - / °?( x) dF x .x -<^b F x -b - 0; 

J x =0 ° ° 


d^T 

dA = kt dA v 
X ° dx 2 X » 


(CIO) 


Substitute the following dimensionless variables: 


X n = X = x ; e = 


And note that 


d 2 T(x) _ T b d 2 g(X) 


dx 


d 2 dX 2 


4 d f L/d 44 4 ktT b d 2 0(X) 

< e eg - J o < e w dF x 0 -x - < F x o -b ■% = y 


oT, 


e 4 (x 


.L/d 


o> ■ f e 4 «> dF X„-X - F X„-b - 


'X=0 


ctT 


bj 


ktT b d 2 6 (X) 


dX 2 
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Thus, 



Equation (Cll) is identical to equation (9), except for the added term q g /oT^, and 
is solved in the same manner to give the temperature along the wall as a function of the 
distance from the opening x/d. The radiant flux q^ emitted from the cell can then be 
determined as in equation (15). 

Performance calculations for the honeycomb absorber at 0. 1 astronomical unit from 
the Sun are compared with collimated solar input conditions in figures 8 and 9 for the 
reference case, which has a base surface temperature T^ of 2060° R (1144 K), a cylin- 
der diameter d of 0. 5 inch (1.27 cm), a cell wall thickness t of 0.003 inch 
(0. 00761 cm), and a thermal conductivity k of 0. 667 Btu per hour-inch-°R (0. 1384 
W/(cm)(K)). As expected, the effect of the added energy input due to the divergent rays 
of the Sun at 0. 1 astronomical unit is to increase the cell wall temperature and, hence, 
the value of e_ relative to the collimated solar input condition. In figure 8, for 

cL 



Figure 8. - Comparison of cell-wall temperature ratio for 
collimated and divergent solar input conditions for 
0. 1-astronomical unit distance from Sun. Reference 
conditions: base temperature, 2060° R (1144 KF, 
diameter, 0.5 inch (1.27 cm); thermal conductivity, 
0.667 Btu per hour-inch-°R (0. 1384 W/(cm)(K)} ; cell 
wall thickness, 0.003 inch (0.00761 cm). 
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Figure 9. - Apparent hemispherical emittance as function of 
cell length -to-diameter ratio. Reference conditions: base 
temperature, 2060° R (1144 K); diameter, 0.5 inch 
(1.27 cm); thermal conductivity, 0.667 Btu per hour- 
inch-°R (0. 1384 W/<cm)(K)) ; cell wall thickness, 0.003 inch 
(0.00761 cm). 


example, at an L/d of 5, the temperature at the cell opening is 0. 535 for the 
collimated case, while for the 0. 1 -astronomical unit case, the temperature increases to 
0. 620 Tjy This increase in wall temperature near the cell opening results in a sub- 
stantial increase in € . For example, in figure 9, at an L/d of 5. 75, e is 0. 1 for 
the collimated case, whereas "e is increased to 0. 19 for the 0. 1 - astronomical unit 
case. 
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